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ABSTRACT 

Material  interfaces  are  common  in  heat  flux  instrumentation  and 
in  modern  thermal  protection  systems.  Interfaces  cause  mis¬ 
matches  which  can  significantly  affect  convective  heat  flux.  A 
literature  survey  on  chemical  and  thermal  mismatches  with  empha¬ 
sis  on  experimental  data  is  presented.  Data  indicate  transient 
variations  in  heat  transfer  coefficient.  Unsteady  boundary  layer 
equations  are  salved  computationally  using  upwind  differences  and 
parabolic  grid  generation.  Results  are  compared  with  data,  and 
other  computational  and  analytical  solutions,  to  investigate 
solution  accuracy  across  a  mismatch.  By  coupling  wall  condi¬ 
tions,  transient  variations  in  heat  transfer  coefficient  at  a 
given  deflection  angle  are  computed  on  Space  Shuttle  Orb iter 
materials  near  an  interface  during  unsteady  pitch  manuevers. 

c„  Specific  heat 

C.,  Skin  friction  coefficient 

h  Heat  transfer  coefficient;  -  T„) 

H  Nondimensional  total  enthalpy 

L  Reference  Length 

M  Mach  number 

p  Nondimensional  pressure 

*  AIAA  Member.  Former  Associate  Professor  of  Aeronautics  & 
Astronautics  at  Air  Force  Institute  of  Technology. 
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Prandtl  number 
Heating  rate 

Nondimensional  surface  radius 
Reynold’s  number 

Stanton  number;  q/  (T0-  Tw)/c^puos 
Nondimensional  time 
Temperature 

Nondimensional  tangential /axial  velocity 
Nondimensional  normal  velocity 
Transformed  normal  velocity 
Nondimensional  tangential /axial  coordinate 
Nondimensional  normal  coordinate 
Transformed  normal  coordinate 
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Transformed  coordinates 
Nondimensional 1 zed  viscosity 
Nondimensional i zed  density 
Deflection  angle 


Subscripts 

<  *  >  Denotes  differentiation 

aw  Denotes  adiabatic  wall 

j  Denotes  normal  grid  point  index 

o  Denotes  stagnation  condition 

T  Denotes  turbulent  quantity 

w  Denotes  wall  quantity 

00  Denotes  freestream  quantity 


Suptmri  Bt 

Two-dimensional /axisymmetric  indicator 
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INTRODUCTION 


Modern  thermal  protection  systems  -for  reentry  vehicles,  such  as 
Space  Shuttle,  usually  have  different  materials  which  interface 
at  numerous  locations.  Heat  flux  instrumentation  imbedded  in 
these  vehicles  and  in  models  used  in  ground  tests  have  similar 
interfaces.  These  interfaces  cause  mismatches  which  may  signifi¬ 
cantly  affect  wall  temperatures  and  convective  heat  flux. 

Several  types  of  mismatches  can  occur  at  material  interfaces, 
as  illustrated  by  the  following  examples.  First,  a  physical 
mismatch  can  be  caused  by  a  forward  or  backward  facing  step  at 
the  surface  o  an  interface.  A  step  causes  flow  disturbances 
which  cause  local  changes  in  convective  heat  flux.  Second,  a 
radiative  mismatch  can  be  caused  by  different  surface  emissivi- 
ties.  Thus,  heat  radiated  from  each  side  of  the  interface  would 
be  different,  and  could  result  in  a  wall  temperature  discontinui¬ 
ty  referred  to  as  a  thermal  mismatch.  A  thermal  mismaten  can 
also  be  caused  by  materials  with  different  thermal  properties, 
cause  different  wall  temperature  responses,  and  result  in  a 
discontinuity  in  the  wall  temperature  distribution  if  conductivi¬ 
ties  are  very  low.  A  fourth  type  of  mismatch,  a  chemical  mis¬ 
match,  can  be  caused  by  coatings  with  different  catalytic  effi¬ 
ciencies  for  the  recombination  of  dissociated  nitrogen  and  oxygen 
(a  chemical  nonequi 1 ibrium  effect).  This  discontinuity  in  sur — 
face  catalysis  can  cause  a  discontinuity  in  recombination  ener — 
gies  which  contribute  to  the  heat  flux.  Of  these  mismatches,  the 
thermal  and  chemical  mismatch  have  often  been  neglected  perhaps 
because  they  were  not  significant,  and  perhaps  because  of  diffi- 
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culty  in  real  i  stical  1  y  accounting  -for  them. 

The  thermal  and  chemical  mismatch  are  difficult  to  account  for 
because  of  dependence  upon  local  and  upstream  wall  conditions. 

For  example,  consider  an  interface  with  a  hot-noncatalytic  wall 
upstream  and  with  a  cold— fully  catalytic  wall  downstream.  During 
reentry,  part  o-f  the  kinetic  energy  of  the  freestream  air  is 
converted  and  stored,  thermally  by  increasing  the  gas  temperature 
and  chemically  by  dissociating  the  gas.  Since  less  energy  is 
transported  to  the  upstream  wall  (heat  -flux)  than  -for  a  coid- 
fuliy  catalytic  wall,  additional  energy  is  convected  downstream 
across  the  interface.  Downstream  of  the  interface,  there  is  a 
jump  in  heat  flux  because  the  wall  is  cold  and  fully  catalytic. 
The  additional  energy  convected  across  the  interface  causes  an 
additional  overshoot  in  heat  flux.  Thus,  the  heat  flux  depends 
on  upstream  and  local  wall  conditions,  which  are  usually  tran¬ 
sient  during  ground  tests  and  reentries. 

Since  numerous  materials  with  different  properties  and  var¬ 
ying  degrees  of  catalytic  efficiency  are  being  used,  a  method  ot 
evaluating  the  significance  of  an  interface  on  the  aerothermody- 
namic  design  of  a  vehicle  (and  heat  flux  instrumentation)  is 
needed.  Experimental  flight  and  ground  tests,  approximate  analy¬ 
tical  solutions,  and  computational  solutions  are  possibilities. 
Only  a  few  references  (which  will  be  presented  in  the  next  sec¬ 
tion)  are  available  on  mismatches  with  one  exception;  numerous 
approximate  analytical  methods  are  published  on  a  thermal  mis¬ 
match  with  a  discontinuous  step  (and  will  not  be  presented).  A 
solution  for  a  chemical  mismatch  with  two  species  is  also  availa- 
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bie,  but  is  not  extended  in  general  for  reacting  mul ti -speci es 
with  finite  recombination  rates.  Some  ground  test  data  are 
available  for  a  thermal  mismatch  with  a  discontinuous  step  on 
simple  conf igurations,  but  ground  facilities  can  not  simulate 
some  variations  in  wall  conditions,  nor  duplicate  flight  condi¬ 
tions.  Computational  solutions  are  becoming  widespread,  and  can 
be  as  general  as  needed,  including  transient  solutions.  But,  the 
accuracy  of  computational  solutions  is  of  concern  because  o+  the 
discontinuities  at  mismatches.  However,  computational  solutions 
have  been  successfully  applied  to  available  ground  and  flight 
data  at  chemical  mismatches  to  correlate  rate  coefficients.  Un¬ 
fortunately,  no  solutions  have  been  computed  to  explain  some 
transient  effects  near  a  mismatch  during  unsteady  pitch  manuevers 
in  bath  ground  and  flight  tests. 

Besides  evaluating  the  significance  of  mismatches,  the  primary 
purpose  of  this  article  is  to  present  and  apply  a  unique  (but 
simple)  computati anal  method  to  investigate  steady  and  unsteady 
solution  accuracy  at.  mismatches  using  a  combination  o+  the  above 
methods.  First,  a  literature  survey  is  presented  which  covers 
available  experimental,  analytical,  and  computational  investiga¬ 
tions  of  the  thermal  and  chemical  mismatch.  Next,  a  computa¬ 
tional  solution  of  the  unsteady  boundary  layer  equations  has  been 
developed  using  three-point  upwind  differences  and  parabolic  grid 
generation.  This  unique  approach  minimizes  computer  resources 
needed  for  transient  solutions,  but  uses  differencing  techniques 
and  arbitrary  grids  applicable  for  Navi er-Stokes  solutions. 

Next,  upwind  difference  solutions  are  validated  by  comparing  with 


analytical  solutions,  central  difference  solutions,  and  experi¬ 
mental  data.  Although  solutions  tor  a  chemical  mismatch  have 
been  computed,  results  are  not  presented  since  a  thermal  mismatch 
is  analogous  and  can  be  validated  easier.  Finally,  the  upwind 
difference  solution  is  coupled  with  wall  temperatures  to  predict 
transient  variations  in  heat  transfer  distributions  during  an 
unsteady  pitch  manuever  with  a  wedge  in  a  wind  tunnel  test. 

LITERATURE  SURVEY 

A  literature  survey  was  conducted  primarily  on  heat  transfer 
data  for  chemical  and  thermal  mismatches,  but  also  included  some 
analytical  and  computational  solutions.  References  for  chemical 
mismatches  will  be  discussed  first.  Then  a  comprehensive  survey 
of  published  experimental  data  for  a  thermal  mismatch  will  be 
presented.  Next,  details  will  be  given  of  some  data  which  will 
be  used  to  compare  with  computational  solutions.  Finally,  some 
data  which  have  indicated  transient  trends  in  heat  transfer 
coefficient  will  be  discussed. 

Chung,  Liu,  and  Mirels  predicted  the  boundary  layer  profiles 
and  heat  flux  for  a  chemical  mismatch  in  a  reacting  binary  mix¬ 
ture.  1  The  overshoot  in  heat  transfer  at  a  chemical  mismatch  was 
confirmed  by  their  theoretical  solution  of  the  boundary  layer 
equations,  although  catalytic  efficiencies  for  various  materials 
were  not  well  known.  Based  on  their  solution,  a  gage  with  known 
catalytic  efficiency  imbedded  in  a  material  ot  unknown  efficiency 
was  proposed  to  determine  catalytic  efficiency.  The  multiple 
species  in  air  during  flight  complicates  the  proposal.  If  air  is 
approximated  as  a  binary  mixture,  then  their  solution  can  vali- 
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date  computational  solutions  across  a  chemical  mismatch. 

Sheldahl  and  Winkler  confirmed  tha  significance  of  a  chemical 
mismatch  in  arc-plasma  tests  by  using  various  combinations  of 
silicon  oxide(SiO)  and  copper (Cu)  coatings  on  a  hemisphere- 
cylinder.*  For  example,  the  heating  on  the  cylinder  with  the  Cu 
coating  was  more  than  twice  the  heating  with  a  SiO  coating.  The 
heating  significantly  quadrupled  with. a  SiO/Cu  interface. 

Stewart,  Rakick,  and  Lanfranco  also  -confirmed  the  significance 
of  a  chemical  mismatch  in  arc-plasma  tests  tor  the  Space  Shuttle 
Qrbiter’s  thermal  protection  system  by  alternatively  using  the 
baseline  coating  of  borosilicate  glass,  and  then  a  catalytic 
overcoat  of  iron-cobalt-chromia-spinal . *  Again,  the  heat  trans¬ 
fer  at  locations  downstream  of  a  noncatalytic/cataiytic  interface 
exceeded  heat  transfer  rates  for  chemical  equilibrium.  Recomoi- 
nation  rate  constants  were  then  adjusted  in  a  computational 
solution  for  a  reacting  boundary  layer  to  match  predicted  ana 
measured  heat  transfer  ratios.  Recombination  rate  constants 
were  then  assumed,  and  similar  overshoots  predicted  for  a  flight 
test  experiment  where  individual  tiles  on  the  lower  surface  of 
the  Qrbiter  were  coated  with  the  same  catalytic  overcoat. 

Shuttle  flight  results  published  by  Rakick,  Stewart,  and  Lan¬ 
franco  (and  by  Scott  and  Derry)  confirmed  the  overshoot  in  heat 
transfer  at  the  chemical  mismatch  caused  by  the  overcoat.*** 
Recombination  rates  had  to  be  adjusted  for  postf light  predic¬ 
tions.  However,  oxygen  and  nitrogen  recombination  rates  can  not 
be  determined  independently  since  their  effects  are  lumped  in 
flight  results.  In  addition,  the  chemical  mismatch  causes  a 
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thermal  mismatch  which  becomes  mo re  significant  at  lower  Mach 
number.  Such  a  combined  chemical /thermal  mismatch  was  assumed  by 
Curry,  Rochelle,  Chao,  and  Ting  for  computational  predictions  at 
the  Orbi tar's  nose  cap/tile  interface.* 

A  comprehensive  survey  of  published  experimental  data  for  a 
thermal  mismatch  provided  several  references.  Data  for  a  thermal 
mismatch  with  a  step  in  wall  temperature  are  given  by  Sees a  for 
incompressible,  turbulent  flow  over  a  flat  plate.7  A  more  thor¬ 
ough  study  for  incompressible,  turbulent  flow  over  a  olate  witn  a 
step  in  wall  temperature  is  given  by  Reynolds,  Kays,  ana  Kline.** 
Laminar  data  for  incompressible  flow  are  given  by  £de  and  Saun¬ 
ders  for  a  flat  plate  in  water,*  but  erratic  data  indicate  subse¬ 
quent  transition  to  turbulent  flow.  McCroskey  published  results 
for  laminar  supersonic  flow  over  a  flat  plate  with  a  step  in  neat 
flux.  »«*•**  Prahara j  and  Foster  presented  flight  data  for  the 
Space  Shuttle  external  tank  during  ascent  which  show  a  signifi¬ 
cant  effect  on  heat  flux  instrumentation  due  to  a  thermal  mis¬ 
match.  **  Data  were  corrected  by  a  factor  of  two.  but  no  correc¬ 
tion  is  available  where  complex  flow  interactions  are  expected. 
Finally,  Durgin  published  excellent  results  for  both  laminar  and 
turbulent  supersonic  flow  over  a  cone. »*• **  Data  from  two  of 
these  references  will  be  discussed  in  more  detail. 

The  experimental  data  tor  turbulent  flow  over  a  fiat  plate  with 
heating  strips  was  selected  to  compare  with  computati onal  solu¬ 
tions.*  The  plate  is  illustrated  m  Figure  1.  The  first  tour 
strips  were  unheated,  and  the  remaining  strips  were  heated  to 
maintain  a  step  in  wall  temperature.  Due  to  the  low  freestream 


B 


Reynold's  number  of  2.86x10*  <0.9x10*  at  th#  thermal  Mismatch),  a 
trip  Mirs  mss  used  after  ths  second  heating  strip  to  transition 
the  flow.  Although  this  data  is  considered  to  be  for  incompress¬ 
ible  flow,  a  freestream  Mach  number  of  0.08  was  computed  for  use 
in  the  compressible  algorithms. 

Mall  temperature  data  on  a  cone  with  a  thermal  mismatch  are 
available  for  compressible  flow  at  Mach  3.5. *■•  **  An 
illustration  of  the  cone  is  shown  in  Figure  2.  The  copper  nose 
was  cooled  for  a  uniform  wall  temperature  of  233" K.  The  aft 
part  of  the  cone  was  made  of  balsa  wood,  which  approximated  an 
adiabatic  wall.  Oata  are  given  for  both  laminar  and  turbulent 
flows  with  freestream  Reynold’s  numbers  of  2.06  and  7.52x10* 
respectively.  Since  numerous  thermocouples  were  along  the  wail 
of  the  cone  (and  uncertainty  avoided  in  determining  heat  flux), 
this  data  was  also  selected  for  computational  comparisons. 

Recent  data  for  a  wedge  shown  in  Figure  3  with  a  thermal  mis¬ 
match  are  also  used  for  computational  comparisons  because  of  some 
interesting  discrepancies. A  steel  wedge  <17.95  x  3d. 46  cm) 
with  inserts  <15.38  x  15.38  cm)  starting  at  17.95  cm  from  tne 
leading  edge  was  pitched  between  3"  to  14"  deflection  angles 
in  a  Mach  14.24  wind  tunnel  to  simulate  Space  Shuttle  uroitmr 
flight  pitch  manuevers.  The  hypersonic  flow  was  laminar  because 
of  low  Reynold's  number  <  l.2xlo*>.  One  insert  was  Fiexidle- 
Reusable-Surf ace- Insulation (FRSI)  from  the  Ore. ter ' •  thermal 
protection  system.  The  heat  transfer  coefficient  to  the  FR8l  was 
dramatically  lower  because  of  the  mismatch,  but  was  also  lower 
than  predicted.  Data  reduction  problems  encountered  by  Moo  indi- 
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cat*  that  th*  h*at  tran*t*r  coefficient  varied  in  time  at  a  given 


deflection  angle  while  th* 


wa*  pitched. **  Thus,  a  hystere¬ 


sis  loop  occurs  when  the  wedge  is  pitched  down  and  then  up.  This 


loop  will  be  predicted  coeputational ly. 

In  another  recent  test  at  liach  6 ,  a  teflon  button  and  coaxial 
thermocouple  gag*  were  imbedded  in  a  blunt-nosed  steel  plat* 
which  was  also  pitched. **  A  hysteresis  loop  was  indicated  by 


ous  data  reduction  methods,  although  results  were  somewhat 


different  from  method  to  method.  Because  of  experimental  diffi¬ 
culties  for  transient  heat  flux  determination,  verification  of 
this  transient  behavior  of  the  heat  transfer  coefficient  by  a 


theoretical  and  computational  solution  is  needed. 


Since  approximate  analytical  methods  can  not  be  extended  in 
general  to  three-dimensional  transient  problems,  computational 


methods  are  becoming 


popular.  The  accuracy  of  heat  transfer 


predictions  are  often  questioned  when  computational  methods  are 


used  especially  for  discontinuities  tsuch  as  a  mismatcm  ana  tor 
coarse  grids.  Fine  grids  can  be  used,  but  computational  re¬ 


source*  may  then  bi 


very  large. 


:ially  if  th*  Navi ar- 


Stokes  equations  tor  th*  flow  are  coupled  with  the  heat  equation 


for  th*  wall  temperatures.  Therefor*,  a  computational  method. 


based  on  upwind  differences,  is  investigated  to  solve  the  un¬ 
steady  boundary  layer  equations  coupled  with  various  boundary 
conditions  for  the  wall  temperatures.  Detail*  of  th*  equations 
and  transf or mat ions  for  arbitrary  grids,  the  upwind  difference 


solution,  and  boundary  conditions  are 


'iced  next. 
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The  unsteady  boundary  layer  aquations  are  solved  computational¬ 
ly  in  terms  of  velocity  and  total  enthalpy.  For  two-dimensional 
flow  <m»0)  or  axisymmetric  flow  <m«l),  the  governing  equations  in 
nondi mensi onal  form  are  given  for  turbulent  flow  as 
R-p%  ♦  <R"*Pu)«  *  <R»Pv>„  -0  < la) 

<pu>*  <pu*)„  *  <puv)v  ■  -p„  ♦  <C  04-*-mt)uv  Jv>/Re  (lb) 

< PH) «  ♦  (puH) H  ♦  (OvH) v  ■  p%  +<C <M/Pr+*T /PrT)Hv  (lc) 

♦  <w (Pr-l ) /P r+nT <PrT-l)/PrT>  <tu»)v]v}/Rs 
where  subscripts  denote  differentiation.  Turbulence  closure  is 
accomplished  with  an  eddy  viscosity  model  following  Cebeci  and 
Smith.  The  total  enthalpy  equation  was  chosen  because  of  a 
truncation  error  criteria.  If  a  Pr andtl  number  of  one  is  as¬ 
sumed.  the  total  enthalpy  is  a  linear  function  of  the  tangential 
velocity  for  constant  properties  and  isothermal  wall  conditions. 
Since  gradients  of  other  variables  (except  density)  are  small 
compared  to  the  tangential  velocity  and  enthalpy  gradients  in  the 
normal  direction,  optimization  of  the  normal  grid  distribution 
for  either  the  velocity  or  total  enthalpy  will  tend  to  minimize 
truncation  error.  Density  gradients  are  usually  not  a  problem  if 
conservative  differencing  or  a  transf or mat  1  on  is  used. 

The  compressible  boundary  layer  equations  are  transformed  to 
eliminate  density  as  a  solution  variable.  Although  usually  used 
in  a  more  general  form,  the  transf ormat  1  on  given  by 

V  -  Jpdy  (2a) 

V  -  v  ♦  V.  ♦  uV„  <2b  > 

is  attributed  to  Dorodnitsyn  by  Kim  and  Chang  who  solved  the 
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unsteady  boundary  layer  equations  using  central  differences. ** 
Some  details  of  the  transf ormation  and  this  approach  are  availa¬ 
ble  in  theses  by  Lange  for  laminar  flow,*0  by  Coffee  for  an 
adaptive  grid  study,**  and  by  Chen  for  turbulent  flow.** 

Contrary  to  most  boundary  layer  solutions  which  use  analytical 
transf oreations,  a  general  transf ormation  for  arbitrary  grids  is 
used  in  addition  to  the  above  transf ormation.  Arbitrary  grids 
are  generated  numerically  by  using  a  new  technique  of  parabolic 
grid  generation,**  by  using  analytical  transf ormation*  for  simi¬ 
lar  solutions,  or  by  using  a  combination  of  techniques.  This 
approach  is  more  flexible,  and  is  also  applicable  to  more  general 
solutions  of  the  Navier — Stakes  equations.  A  general  transforma¬ 
tion  for  arbitrary  grids  in  x-Y  space  to  space  (corresponding 

to  a  square  grid  of  unit  step  size  in  a  rectangular  domain)  is 
then  applied  to  Eqs.  <1>  using  the  chain  rule.  One  metnc(  >:n>, 
corresponding  to  normal  coordinate  lines,  must  be  zero  to  be 
consistent  with  traditional  boundary  layer  assumptions.  The 
resulting  equations  are  given  by 

(R-V)  -  <Y£<R-u)n  -  Yp  (R**u>  ^1/x^  (3a) 

u%  ♦  uu./K>  <V  -  Y_u/x_)u  /Y  ■  -p_/x_  (3b) 

*  *  t  t  n  n  t  (H+uT)uTi/Yn  ln/ (YnRe) 

W%  ♦  uH^ / x ^  ♦  (V  -  *£u/>;£  >Hp/Yp  •  tp  t»/Pr+*r/Pr  T)Hp /Yp  3n  /  ( YpRe) 

♦  P%/p  ♦  CpC*<(Pr-i  > /Pr^HT(PrT-l) /PrTJ  (tu*)p/Yp>p/ (YpRe) 

(  3c ) 

in  a  form  convenient  for  computational  solution. 

The  upwind  difference  solution  which  was  developed  is  different 
from  usual  approaches  pr i atari  ly  in  three  ways.  First,  upwind 
differences  using  either  a  three-point  forward  or  backward  dif- 
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ference  (depending  on  the  signs  of  Appropriate  eigenvalues  of 
diagonal  matrices)  are  used.  This  was  suggested  for  the  boundary 
layer  equations, **  but  is  analogous  to  the  approach  used  by 
Hodge,  Stone,  and  Miller  previously  in  the  solution  of  the  incom¬ 
pressible  Navi er — Stokes  equations. *•  Second,  the  metrics  of  the 
transf ormation  are  calculated  in  numerous  ways.  Third,  the  dif¬ 
ference  solution  is  simplified  because  of  the  upwind  differences. 

Upwind  differences,  instead  of  traditional  central  differences, 
are  used  for  the  nonlinear  convective  terms  in  Eqs.  (3b,  3c ) 
because  of  several  advantages.  Linearized  matrix  equations  are 
more  diagonally  dominant  when  upwind  differences  are  used,  and 
are  unconditionally  stable  for  Eqs.  (3b,  3c).  Artificial  viscos¬ 
ity  is  not  necessary  since  the  truncation  error  is  diffusive,  as 
apposed  to  dispersive  for  central  differences.  The  magnitude  of 
the  truncation  error  term  for  an  upwind  difference  is  larger  by  a 
factor  of  two,  but  is  second  order  also  as  shown  below. 

Truncation  error  depends  on  the  grid  distribution,  and  how  the 
metrics  are  calculated.  If  the  grid  is  generated  numerically, 
then  the  metrics  are  approximated  by  finite  differences  causing  a 


>tr  ic 

error.  For 

ex ample 

if  V  O, 

then 

VUy 

m 

MVNCfc (3Uj 

—  4u,.| 

♦  u4*»> 

♦  T.E. 3/ 

-  y4_»> 

♦  M.E. 3 

<  4a  > 

M.E. 

m 

“vnnn/4  ♦ 

•  •  a 

(4b  > 

T.E. 

- 

unnn/3  ♦ 

s  •  a 

(4c) 

where  M.E.  refers  to  "metric  error"  and  T.E.  to  "truncation 
error".  The  grid  distribution  could  be  optimized  so  that  T.E.  is 
zero,  but  M.E.  could  be  large.  Another  option  would  be  to  gene¬ 
rate  the  grid  analytically  so  the  metrics  could  be  calculated 
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exactly,  but  T.E.  could  be  large  for  a  coarse  grid.  A  combina¬ 
tion  of  options  is  used  however,  and  is  very  useful  near  singu¬ 
larities.  For  example,  the  metric  is  calculated  analytically 
according  to  similar  solutions  near  a  leading  edge  singular  ity, 
Mhile  downstream  of  the  singularity  the  grid  may  be  controlled 
numerically  to  minimize  truncation  error.  Thus,  metric  calcula¬ 
tion  is  flexible,  and  independent  of  the  difference  solution. 

Solution  of  the  upwind  difference  equations  is  simplified  com¬ 
pared  to  central  difference  methods.  Although  matrix  bandwidth 
increases,  iterative  convergence  necessary  for  nonlinear  equa¬ 
tions  is  enhanced  by  upwind  differences  because  of  diagonal 
dominance.  Dominance  is  obtained  even  for  flows  where  convec¬ 
tion  dominates,  since  the  convective  terms  in  Eqs.  (3b,  3c >  are 
approximated  by  one-sided  backward  or  forward  differences  depend¬ 
ing  on  the  sign  of  their  coefficients  which  depend  on  the  pre¬ 
vious  iteration.  Dominance  is  further  enhanced  by  a  two-point  or 
three-point  backward  difference  for  time  derivatives.  Because 
of  diagonal  dominance,  a  Successive-Over — Rel ax at ion (SOR)  or  tri- 
di agonal  matrix  solution  technique  for  Eqs.  (3b,  3c)  is  suf¬ 
ficient.  However,  a  quad-diagonal  matrix  solution  technique  is 
u*******  Eqs.  <2,  3a)  are  then  integrated  by  trapezoidal  rule.  A 
few  iterations  are  necessary  for  steady  or  unsteady  solutions, 
and  converge  rapidly  if  wall  boundary  conditions  are  simple. 

boundary  conditions  are  specified  or  coupled  with  the  above 
finite  difference  solution.  Because  of  the  parabolic  nature  of 
the  boundary  layer  equations,  the  normal  velocity  is  specified  at 
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the  wall  and  computed  at  the  outer  boundary.  The  physical  coord¬ 
inate  (y)  is  likewise  specified  at  th®  wall  and  than  computed  at 
th®  outer  boundary  (for  external  problems).  An  inv^scid  solution 
provides  the  outer  edge  conditions  for  the  velocity,  enthalpy, 
and  pressure.  The  pressure  across  the  boundary  layer  is  known 
since  the  normal  pressure  gradient  is  zero  from  the  normal  mo- 
eentum  equation.  The  no  slip  condition  is  prescribed  at  the 
wall,  and  only  the  wall  enthalpy  is  unknown. 

Various  conditions  are  used  for  wall  enthalpy.  One  option  is 
to  specify  the  wall  temperature.  An  isothermal  wall,  or  a  wall 
temperature  step  corresponding  to  the  experimental  data  on  the 
plate*  and  wedge,1*  are  obvious  choices.  Specified  wall  tempera¬ 
tures  fallowed  by  an  adiabatic  wall  is  another  option  correspond¬ 
ing  to  experimental  data  on  the  cone.1'*  A  radiation  equilibrium 
condition  with  a  resistance-capacitance  analogy  to  simulate  con¬ 
duction  effects  is  used  for  the  insert  in  the  wedge. ao  wlqo- 
rithms  for  the  one-di men si onal ,  two-dimensional,  or  axisymmetric 
heat  equation  for  the  wall  can  be  coupled  to  the  solution,  and 
heat  fl  ux  gages  eventually  investigated.  However,  the  boundary 
conditions  were  simplified  to  investigate  the  accuracy  of  the 
computational  solutions  and  significance  of  mismatches. 

RESULTS 

The  accuracy  of  the  upwind  difference  algorithm  is  validated  by 
comparisons  with  analytical  solutions  and  other  algorithms  using 
central  differences.  The  accuracy  of  the  difference  solutions 
across  a  thermal  mismatch  and  the  significance  of  mismatches  are 
then  implied  by  additional  comparisons  with  experimental  data. 


IS 


Finally,  unique  results  for  an  unsteady  problem  are  presented. 
Results  are  computed  for  three  simple  geometric  shapes  to  avoid 
additional  complexities  caused  by  curvature.  These  shapes  in¬ 
clude  a  flat  plate,  a  cone,  and  a  wedge. 

Flat  Plate 

First,  the  skin  friction  Mas  computed  on  a  flat  plate  using  the 
upwind  difference  algorithm,  since  an  exact  solution  is  available 
for  incompressible  flow.  Since  the  algorithms  are  for  compressi¬ 
ble  flow,  a  Mach  number  of  0.001  and  constant  temperature  were 
assumed.  The  same  30x30  grid  generated  by  a  central  difference 
algorithm  with  an  analytical  transf ormat i on  was  used.  The  met¬ 
rics  were  very  sensitive  to  how  they  were  approximated,  and 
therefore  the  metric  was  calculated  anal yti cal  1 y .  The  skin 
friction  was  less  than  the  exact  theoretical  solution  by  no  more 
than  one  percent,  whereas  the  central  difference  algorithm  had  no 
•ore  than  a  one-half  percent  overshoot.  Thus,  the  accuracy  tor 
incompressible  flow  is  validated. 

Other  grids  were  also  used  since  difference  algorithms  are  very 
grid  dependent.  The  upwind  algorithm  did  not  diverge  when  a  very 
coarse  grid  was  used,  but  the  central  algorithm  did.  Both  algo¬ 
rithms  demonstrated  second-order  convergence  when  the  step  size 
decreased.  Second  order  upwind  differences  were  usually  neces¬ 
sary,  however,  since  the  skin  friction  with  first  order  upwind 
differences  for  the  above  grid  was  only  within  eight  percent. 
Similar  results  were  also  obtained  for  a  thermal  mismatch. 

Next,  the  upwind  algorithm  was  validated  for  turbulent,  incom¬ 
pressible  flow  over  a  plate  with  a  thermal  mismatch.  The  heat 


transfer  for  the  unheated/heated  plats  in  Figure  1  with  a  11“K 
temperature  step  was  computed,  and  agrees  well  with  experimental 
data  as  shown  in  Figure  4.  Predictions  -for  laminar  and  turbulent 
flow  with  and  without  the  mismatch  are  also  shown.  Note  that  the 
percentage  change  in  heat  transfer  due  to  the  thermal  mismatch  is 
greater  for  laminar  flow.  The  accuracy  was  more  sensitive  to 
turbulence  modeling  than  to  the  numerics  for  a  61x59  grid  with 
several  grid  points  in  the  sublayer.  Therefore,  computational 
accuracy  for  a  mismatch  is  addressed  better  using  laminar  data. 
CSDS 

Since  steady  wall  temperature  data  are  available  tor  both 
laminar  and  turbulent  flow  over  a  cone  with  a  thermal  mismatch,  1  ■* 
upwind  difference  solutions  were  computed  to  validate  accuracy 
across  a  mismatch.  If  the  coaled  nose  of  the  cone  is  assumed  to 
have  an  isothermal  wall  and  no  boundary  layer /flow  interaction, 
then  the  heat  transfer  prediction  is  analogous  to  a  flat  plate 
and  wedge  according  to  boundary  layer  theory,  except  for  an 
increase  in  heat  transfer  by  a  factor  of  /2.  Thus,  the  computed 
heat  transfer  coefficients  at  the  nose  can  be  compared  with 
analytical  solutions.  Then,  the  computed  adiabatic  wall  tempera¬ 
ture  distribution  downstream  of  the  mismatch  is  compared  to  the 
measured  wall  temperature  distribution. 

The  isothermal  heat  transfer  coefficient  on  the  cool ed  nose, 
predicted  by  the  upwind  difference  solution  for  compressible  flow 
with  a  fine  grid  (61x59),  compared  with  the  central  difference 
solution  to  the  third  decimal.  The  magnitude  of  the  heat  trans¬ 
fer  coefficient  also  compared  within  a  few  percent  to  approximate 
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analytical  solutions  depending  on  the  viscosity  law.  Further 
comparisons  are  not  necessary  since  Reynold’s  analogy  leads  back 
to  the  skin  friction  comparisons  above  for  incompressible  flow 
over  a  flat  plate,  it  the  viscosity  is  assumed  constant. 

Next,  adiabatic  wall  temperatures  downstream  of  the  thermal 
mismatch  were  computed  and  compared.  Durgin  obtained  excellent 
comparisons  between  the  data  and  approximate  analytical  solu¬ 
tions,  14  if  actual  wall  temperatures  were  used  upstream  o+  the 
mismatch.  The  measured  and  predicted  wall  temperature  distribu¬ 
tions  along  the  surface  of  the  cone  are  shown  in  Figure  5  tor 
both  laminar  and  turbulent  flow.  Good  agreement  is  shown  consid¬ 
ering  that  a  uniform  temperature  was  assumed  on  the  cooled  nose. 
In  fact,  this  illustrates  that  the  error  is  no  greater  than  the 
experimental  error  in  maintaining  an  isothermal  wall  on  the  nose. 

If  the  analytical  and  computati anal  solutions  are  amplified 
just  downstream  of  the  thermal  mismatch  and  discontinuity,  how¬ 
ever,  a  discrepancy  is  evident.  The  jump  in  heat  transfer  or 
wall  temperature  at  the  mismatch  is  smeared  by  numerical  trunca¬ 
tion  error,  and  is  smeared  by  conduction  effects  in  any  experi¬ 
ment.  Therefore,  the  required  grid  size. at  the  mismatch  depends 
on  the  material  properties.  To  account  for  these  properties,  and 
to  also  compute  wall  temperatures  in  general,  the  computational 
solution  may  need  to  be  coupled  with  a  computational  solution  of 
the  heat  equation  for  the  materials.  Since  the  wail  temperatures 
are  usually  transient,  the  unsteady  equations  are  solved. 

wsdgs 

Several  algorithms  were  used  to  predict  the  heat  transfer  on 
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the  wedge  shown  in  Figure  3  with  a  thermal  mismatch.  Approximate 
analytical  solutions,  central  difference  solutions  of  the  steady 
boundary  layer  equations,  solutions  of  the  Navi er-Stokes  equa¬ 
tions  using  an  explicit  Mac Cor mack ' s  algorithm,  upwind  difference 
solutions  of  the  steady  boundary  layer  equations,  and  upwind 
difference  solutions  of  the  unsteady  boundary  1 ayer  equations 
have  been  compared  to  limited  exper 1  mental  data  on  the  wedge  tor 
constant  deflection  angle  and  for  a  pitch  manuever.  Steady 
results  for  a  thermal  mismatch  corresponding  to  this  data  have 
been  obtained  by  assuming  an  appropriate  wall  temperature  step. 
The  analytical  solutions  and  steady  boundary  layer  results  are 
given  first.  Since  the  flow  was  hypersonic  and  some  viscous- 
shock  interaction  was  expected,  solutions  of  the  Navier-Stokes 
equations  are  then  presented.  Finally,  since  transient  effects 
have  been  postulated,  upwind  difference  solutions  of  the  unsteady 
boundary  1 ayer  equations  coupled  with  a  wall  condition  for  a 
thermal  mismatch  are  presented  for  a  pitching  wedge. 

Several  steady  solutions  of  the  boundary  layer  equations  for 
the  wedge  with  a  specified  wall  temperature  step  of  95“ K  were 
computed.  At  14"  deflection  angle,  the  ratio  of  the  heat 
transfer  coefficient  to  a  reference  coefficient  at  0“  is  shown 
in  Figure  6.  Computati onal  and  approximate  analytical  solutions 
compare  as  mentioned  above  for  laminar  flow  on  a  cone.  Although 
predictions  account  for  most  of  the  significant  decrease  in  heat 
transfer  coefficient.,  the  experimental  heat  transfer  is  lower 
near  the  thermal  mismatch,  but  recovers  and  agrees  downstream. 

Since  a  strong  viscous-shock  interaction  is  possible  at  hyper- 
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sonic  speeds  over  the  -front  one— third  o-f  the  wedge  (especially  at 
a  lower  deflection  angle  of  three  degrees),  the  two-dimensional 
Navier->Stokes  equations  were  solved  computationally  by  Roberts 
using  a  62x30  exponentially  stretched  grid.a<*  Results  for  both 
an  isothermal  wall  and  a  thermal  mi smatch  with  an  assumed  temper¬ 
ature  step  of  88“  K  are  shown  in  Figures  7  and  8.  The  ratio  of 
the  wall  pressure  to  freestream  pressure  is  compared  in  Figure  8 
to  the  inviscid  solution,  which  is  used  by  boundary  layer  theory. 
The  wall  pressure  is  higher,  as  expected  because  of  'viscous-shock 
interaction.  The  thickening  of  the  thermal  layer  downstream  of 
the  thermal  mismatch  then  caused  another  interaction,  which  is 
shown  by  a  small  increase  in  wall  pressure  further  downstream. 
However,  the  interactions  only  cause  small  increases  in  the  heat 
transfer  as  shown  in  Figure  8.  The  heat  transfer  prediction  is 
also  higher  than  the  experimental  data. 

Since  time  variations  in  heat  transfer  coefficient  have  been 
postulated,  unsteady  effects  and  solutions  were  investigated. 
Numerous  effects,  although  perhaps  small,  cause  time  variations 
in  the  wind  tunnel.  Among  these  are  variations  in  tunnel  condi¬ 
tions,  increasing  wall  temperatures  near  the  leading  edge,  con¬ 
duction,  and  radiation.  An  advantage  of  computational  methods  is 
that  these  effects  can  be  investigated  separately.  Therefore, 
the  radiation  equilibrium  boundary  condition  with  the  resistance- 
capacitance  analogy  was  used  downstream  of  the  stainless  steel/ 
FRSI  interface  to  simulate  thermal  lag  due  to  conduction.  Since 
the  flow  is  hypersonic,  a  quasi-steady  approach  (steady  inviscid 
solutions)  is  used  for  edge  conditions,  and  pitch  rate  effects 
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neglected.  The  unsteady  boundary  layer  equations  were  then  cou¬ 
pled  with  wall  boundary  conditions,  and  solved  tor  a  simulated 
injection  and  pitch  manuever  of  the  wedge  in  the  wind  tunnel. 

The  wall  temperature  distribution  was  thus  obtained  downstream 
of  the  thermal  mismatch  at  each  time  step.  After  injection  tc  a 
constant  14"  deflection  angle,  the  wall  temperatures  were  al¬ 
lowed  to  reach  equilibrium.  The  wall  temperature  distribution  at 
equilibrium  is  shown  in  Figure  9  correspond! ng  to  the  highest 
temperatures.  These  temperatures  also  correspond  to  the  start  of 
the  pitch  manuever.  After  the  pitch  down  to  3°,  the  wall 
temperatures  decrease,  but  lag  the  lowest  wall  temperatures  shown 
which  occur  after  the  five  second  hold  at  3“  (and  before  the 
pitch  up  to  14°).  At  the  end  of  the  pitch  manuever,  the  wall 
temperatures  are  less  than  the  equilibrium  temperatures  due  to 
the  lag.  Equilibrium  is  eventually  reached  again  with  the  radia¬ 
tive  heat  rate  equal  to  the  convective  heat  rate. 

The  convective  heat  rate  distribution  was  also  computed  at  each 
time  step,  but  the  postulated  hysteresis  loop  is  best  illustrated 
at  a  given  location  as  shown  in  Figure  10.  The  heat  transfer 
coefficients  at  14*  (and  start  of  the  pitch  manuever)  decrease 
as  the  wedge  is  pitched  down.  During  the  five  second  hold  at 
3°,  the  heat  transfer  increases.  Then  as  the  wedge  is  pitched 
up,  the  heat  transfer  is  higher  than  during  the  pitch  down,  and 
increases  to  a  maximum  at  the  end  of  the  pitch,  manuever.  As  the 
wall  temperatures  increase  towards  equilibrium,  the  heat  transfer 
eventually  decreases  back  to  equilibrium.  It  should  be  noted 
that  this  significant  time  variation  in  heat  transfer  occurs  just 
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downstream  of  a  thermal  mismatch.  Further  downstream,  the  varia¬ 
tion  and  hysteresis  loop  are  smaller.  However,  unsteady  flow 
effects  were  very  small. 

Although  the  external  flow  was  assumed  to  be  quasi -steady, 
unsteady  terms  were  still  included  in  the  boundary  layer.  Some 
additional  lag  in  the  wall  temperatures  and  an  overshoot  in  heat 
transfer  could  be  caused  due  to  finite  convection  velocities  near 
the  wall.  Solutions  with  and  without  the  unsteady  terms  were 
computed,  and  temporal  truncation  error  was  checked  by  using 
smaller  time  steps.  A  very  small  overshoot  in  the  heat  transfer- 
ratio  of  less  than  one-tenth  percent  was  obtained  when  a  ,er , 
small  time  step  of  25  milliseconds  was  used.  buch  a  small  ertgc: 
is  expected  for  hypersonic  flow  where  only  a  very  small  fraction 
of  the  boundary  layer  is  subsonic.  The  time  scale  for  the  hyper¬ 
sonic  "flow  is  essentially  three  or  more  orders  of  magnitude  less 
than  the  time  scale  for  heat  conduction  at  the  wall.  Thus,  a 
completely  quasi -steady  approach  using  steady  boundary  layer 
solutions  or  approximate  analytical  solutions  would  be  accurate. 
For  lower  velocities  and  higher  pitch  rates,  unst’eady  solutions 
may  be  required. 

CONCLUSIONS  AND  RECOfUENDAT I QNS 

The  flexibility  of  algorithms  using  upwind  differences  and  a 
general  transf ormation  has  been  demonstrated  due  to  improved 
diagonal  dominance.  Steady,  unsteady,  and  quasi -steady  solutions 
were  computed  efficiently,  and  wall  boundary  conditions  with 
dramatically  different  time  scales  were  coupled  with  the  solu¬ 
tion.  Although  alternative  solutions  are  available  tor  a  boundary 


layer,  the  unsteady  nature  of  viscous  solutions  of  the  Navier- 
Stokes  equations  may  necessitate  a  similar  approach.  This  flexi¬ 
ble  approach  is  also  applicable  to  Navi er-Stokes  solutions. 

The  accuracy  of  an  upwind  algorithm  for  the  unsteady  boundar y 
layer  equations  was  also  validated  by  comparisons  with  analytical 
solutions,  central  difference  solutions,  and  #;.per  i  mental  data. 
The  error  in  skin  friction  and  heat  transfer  was  shown  to  be 
larger  for  upwind  differences  than  for  central  differences,  fc"t 
of  the  same  order  magnitude.  Explicit  artificial  viscosity, 
which  often  degrades  the  accuracy  of  central  difference  solu¬ 
tions,  was  not  neccessary  with  upwind  differences  even  for  verv 
coarse  grids.  However,  second-order  upwind  differences  were 
necessary  with  coarse  grids,  especially  for  mismatches. 

The  accuracy  of  difference  methods  across  a  thermal  mismatch, 
and  significance  of  mismatches,  was  also  demonstrated.  Second  - 
order  solutions,  but  not  first-order,  were  snown  to  be  accurate 
across  a  mismatch.  The  solutions  ana  tne  experi  mental  data  also 
demonstrated  that  the  heat  transfer  can  be  si gm 1 1  cant i y  attested 
locally  by  mismatches.  The  dependence  on  thermal  properties  and 
wall  boundary  conditions  were  demonstrated  by  predicting  time 
variations  in  the  heat  transfer  dunno  a  pitch  manuever  6t  a 
wedge  with  a  thermal  mismatch. 

Since  modern  vehicles  and  heat  flux  i nstrumentat 1  on  have  numer¬ 
ous  material  interfaces  where  thermal  and  chemical  mismatenes  are 
possible,  some  method  should  be  used  to  evaluate  the  significance 
of  the  mismatch.  Analytical  approx l mat 1 ons  can  often  be  used, 
and  appropriate  wall  boundary  conditions  should  be  coupled  with 
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